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ABSTRACT 

We report the discovery of a 206 ms pulsar associated with the TeV y-ray source HESS J1 640—465 using the 
Nuclear Spectroscopic Telescope Array ( NuSTAR ) X-ray observatory. PSR J1 640— 4631 lies within the shell- 
type supernova remnant (SNR) G338.3— 0.0, and coincides with an X-ray point source and putative pulsar wind 
nebula (PWN) previously identified in XMM-Newton and Chandra images. It is spinning down rapidly with period 
derivative P = 9.758(44) x 10 -13 , yielding a spin-down luminosity E = 4.4 x 10 36 erg s _1 , characteristic age 
z c = P/2P = 3350 yr, and surface dipole magnetic field strength B s = 1.4 x 10 13 G. For the measured distance of 
12 kpc to G338.3— 0.0, the 0.2-10 TeV luminosity of HESS J1640— 465 is 6% of the pulsar’s present E. The Fermi 
source 1FHL J1640.5— 4634 is marginally coincident with PSR J1640— 4631, but we find no y-ray pulsations in a 
search using five years of Fermi Large Area Telescope (LAT) data. The pulsar energetics support an evolutionary 
PWN model for the broadband spectrum of HESS J1640— 465, provided that the pulsar’s braking index is n ~ 2, 
and that its initial spin period was Pq ~ 15 ms. 

Key words: ISM: individual objects (G338.3-0.0, HESS J1640-465, 1FHL J1640.5-4634, AX J1640.7-4632, 
XMMU J164045.4— 463131) - ISM: supernova remnants - pulsars: individual (PSR J1640-4631) - stars: neutron 
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1. INTRODUCTION 

The detection by the HESS Galactic plane survey (Aharonian 
et al. 2005) of 10 12 eV radiation coming from the diffuse remains 
of supernovae has opened up a new window on the evolution 
of these energetic stellar remnants. More than two-thirds of the 
>80 Galactic TeV sources 12 are supernova remnants (SNRs) or 
pulsar wind nebulae (PWNe), the latter being the most numerous 
class. 

High-energy radiation from PWNe is produced as a result of 
the interaction of the pulsar wind with the surrounding medium. 
Many young PWNe are found inside shell SNRs, where the 
emission begins at a termination shock situated close to the 
pulsar. Downstream of the shock the relativistic electrons radiate 
synchrotron photons from radio through X-rays. The same 
electrons up-scatter ambient photons, creating a second, broad 
spectral peak in the y- ray band. Some SNR shells also emit 
high-energy y- rays, although the dominant mechanism there 
is less clear. One possibility is the hadronic scenario, decay 
of 7T° mesons created in collisions between shock accelerated 
cosmic-ray protons and thermal gas. However, leptonic models 
involving inverse Compton scattering of ambient photons are 
also invoked (for a review, see Reynolds 2008). 

HESS J1640— 465 (Aharonian et al. 2006) is coincident 
with G338.3— 0.0, a shell-type, 8' diameter SNR (Shaver & 


11 Affiliate Member. 

12 TeVCat, http://tevcat.uchicago.edu/. 


Goss 1970; Whiteoak & Green 1996). The TeV emission was 
described as resolved, but centrally peaked. Deeper observations 
(HESS Collaboration 2014) show a more extended source 
than first reported, with Gaussian a = 4 .'3 =b 0'2 and a flux 
above 200 GeV of 1.65 x 10 -11 erg cm -2 s _1 . Based on 
21 cm Hi absorption spectra to G338.3— 0.0 and adjacent Hu 
regions, Lemiere et al. (2009) concluded that the distance to 
G338.3— 0.0 is in the range 12-13.5 kpc, which makes HESS 
J1 640—465 the most luminous TeV source in the Galaxy, with 
L(0.2-10 TeV) = 2.8 x 10 35 (d/12 kpc) 2 erg s“*. 

An XMM-Newton observation by Funk et al. (2007) identi- 
fied a highly absorbed X-ray point source coincident with the 
HESS source with clear indication of extended emission. The 
X-ray components were subsequently well resolved by Chandra 
(Lemiere et al. 2009). No X-ray emission is detected from the 
SNR shell, probably because of its low temperature and large 
intervening column density, Nu ~ 1.4 x 10 23 cm -2 . The power- 
law spectrum of the Chandra point source, with photon index 
r PSR ~ 1.1, is consistent with a pulsar origin. The compact 
X-ray PWN is ^1 '2 in diameter, smaller than the TeV source, 
while its spectral index, Tpwn ~ 2.5, is much steeper than that of 
the pulsar. Its spectrum steepens further in the outer parts, which 
is evidence of synchrotron aging. No radio counterpart to the 
X-ray PWN is detected in high-resolution images of 
G338.3— 0.0 using the GMRT at 235, 610, or 1280 MHz, and 
ATCA maps at 1290 and 2300 MHz, and no radio pulsations 
were found by the GMRT at 610 and 1280 MHz (Castelletti 
et al. 2011). 
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Table 1 

Log of NuSTAR Observations and Period Measurements 


ObsID 

Start Date 
(UT) 

Exposure/ 
Span (ks) 

Rate a 

(s- 1 ) 

Start Epoch 
(MJD) 

Period b 

(s) 

Z 2 

40014016001/17001 

2013 Jun 22 

48.6/96.5 

0.011 

56,465.91194958 

0.206443040(33) 

59 

30002021002/21003 

2013 Sep 29 

89.9/166.4 

0.011 

56,564.29199072 

0.206451335(17) 

79 


Notes. 

a Background subtracted 3-25 keV count rate in a 30" radius aperture. 

b Period derived from a Z\ test. The Monte Carlo derived 1 o uncertainty on the last digits is in parentheses. 


The Fermi source 1FGL J1 640. 8— 4634, coincident in po- 
sition with HESS J1640— 465, was interpreted by Slane et al. 
(2010) as leptonic emission from a PWN. Their derived GeV 
spectrum is not cut off as would be the case for magneto- 
spheric pulsar emission, but is continuous with the HESS spec- 
trum. In contrast, HESS Collaboration (2014) interpreted the 
joint Fermi/ HESS spectrum and revised TeV extent in terms of 
hadronic emission from a portion of the SNR shell interacting 
with dense interstellar gas. More recent analyses of additional 
Fermi data (Nolan et al. 2012; Acero et al. 2013; Ackermann 
et al. 2013) show less flux than Slane et al. (2010) found, and 
the revision will affect both models. 

In this paper, we present the discovery of X-ray pulsations 
from HESS J1640— 465/G338.3— 0.0 in data acquired as part of 
the NuSTAR survey of the Norma Arm region of the Galactic 
plane (Harrison et al. 2013; Fornasini et al. 2014), and a 
follow-up observation to determine the pulsar’s energetics. 
Section 2 describes the NuSTAR observations, pulsar analysis, 
and spectrum. In Section 3 we report on an unsuccessful search 
of the Fermi LAT data for y- ray pulsations. In Section 4 we 
discuss the properties of HESS J1640— 465/G338.3— 0.0 in 
the context of previous leptonic PWN models for its broad- 
band spectrum, and use the measured spin parameters of PSR 
J1 640—463 1 in a revised evolutionary PWN model. Some of the 
implications of our model, and a comparison of its assumptions 
and results with the hadronic model of HESS Collaboration 
(2014), are presented in Section 5. 

2. NuSTAR OBSERVATIONS 

The field containing HESS J1 640—465 was observed by 
NuSTAR in three offset pointings of ~25 ks each on 2013 June 
20-24 as part of a survey of the Norma region of the Galactic 
plane. However, the source fell in the chip gap for the first 
observation (ObsID: 40014011002) so we exclude this data set 
in the final timing and spectral analysis. A follow-up, dedicated 
90 ks observation of HESS J1640— 465 was performed on 2013 
September 29. An observation log is presented in Table 1. 
Data were collected using NuSTAR' s two co-aligned X-ray 
telescopes, with corresponding focal plane modules FPMA 
and FPMB. These telescopes provide 18" FWHM imaging 
resolution over the 3-79 keV X-ray band, with a characteristic 
spectral resolution of 400 eV FWHM at 10 keV (Harrison et al. 
2013). The reconstructed NuSTAR coordinates are accurate to 
7. "5 at 90% confidence. The relative timing accuracy of NuSTAR 
is limited to ~2 ms rms, after correcting for thermal drift of the 
on-board clock, with the absolute timescale shown to be better 
than <3 ms (Mori et al. 2014). 

Data were processed and analyzed using NuSTARDAS vl.3.1 
and HEAS0FT v6.15.1 and the Calibration Database (CALDB) 
files from 2013 August 30. Our analysis used the standard 
level- 1 filtered event files generated by nupipeline. The 



Figure 1. NuSTAR 3-79 keV exposure-corrected image of the field containing 
HESS J1 640— 465 acquired on 2013 June 22. The newly detected pulsar is co- 
incident with CXOU J164043. 5— 463135 (tick marks). Significant 0.5-10 keV 
emission from the PWN surrounding the pulsar is outlined by the Chandra 
contours (red) and enclosed by the spectral extraction region (dotted ellipse). 
SNR G338.3— 0.0 is indicated by the MOST 843 MHz contours (blue; Whiteoak 
& Green 1996). The TeV extent of HESS J1640— 465 is denoted by the solid 
circle (HESS Collaboration 2014). PSR J1640— 4631 is just outside the 95% 
confidence error ellipse of the >10 GeV Fermi source 1FHL J1640.5— 4634 
(dashed circle; Ackermann et al. 2013). 

(A color version of this figure is available in the online journal.) 


observations were free of significant time variable particle 
background contamination. However, each data set has a unique 
spatial background pattern across the field-of-view, which 
affects the sensitivity to a pulsed signal. 

Figure 1 displays the NuSTAR 3-79 keV image of HESS 
J 1640— 465, obtained on 2013 June 22 with the source on- 
axis (ObsID: 40014016001). The image has been corrected for 
exposure, smoothed using a o — 7. "4 Gaussian kernel, and 
scaled linearly. There is one significant source whose position 
is 1" from the putative pulsar CXOU J164043. 5— 463135, 
which is coincident within the NuSTAR positional uncertainty. 
There is evidently diffuse emission surrounding the point 
source, consistent with that found in XMM-Newton and Chandra 
images, and interpreted as a PWN (Funk et al. 2007; Lemiere 
et al. 2009). This, plus its X-ray flux and spectrum (see below) 
leave no doubt that the NuSTAR source is the counterpart of 
CXOU J1 64043. 5— 463135. The most recent Chandra image 
obtained on 2011 June 6 (Fornasini et al. 2014) shows that this 
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Figure 2. NuSTAR discovery of PSR J1640— 4631. Left: superimposed power spectra for the 2013 June 22 discovery observation and the 2013 September 29 follow-up 
showing the significant 4.84 Hz signal and its change in frequency. Right: the sum of background-subtracted pulse profiles from both epochs, each folded using the 
timing parameters in Table 2, and artificially aligned in phase. Phase zero is arbitrary and two cycles are shown for clarity. Data in both panels include 3-25 keV 
photons extracted from a 30" radius aperture. 


source has not varied compared to the two previous Chandra 
observations of 2010 June 19 and 2007 May 11. 

2.7. Timing Analysis 

To search for pulsations from CXOU J 164043. 5— 463135, 
we examined NuSTAR data obtained during the two Norma 
Arm survey pointings, adjacent in time. The day-long data set 
yields a total of 913 photons from the two focal plane modules, 
extracted in the optimal 3-25 keV energy bandpass using a 30" 
radius aperture centered on the source location. Photon arrival 
times were corrected for clock drift and converted to barycentric 
dynamical time (TDB) using the JPL DE200 ephemeris and 
the Chandra coordinates of the point source. The arrival times 
were binned at 2 ms and searched for coherent pulsations up 
to the Nyquist frequency using a 2 26 bin fast Fourier transform 
(FFT). We found a 206 ms coherent signal with a power of 
Pfft = 58, significant at the 99.997% (~4 a) confidence level. 
This motivated a second-epoch observation from which we were 
able to confirm the detection (Pfft = 77) and measure the spin- 
down rate of the pulsar as described next. 

We computed the spin-down rate from the difference in 
frequency over the 100 day time span between the June and 
September observations. For each data set, we generated a 
refined frequency measurement by oversampling the signal 
using the Z 2 test statistic around the known frequency (Table 1). 
We then re-fitted for frequency by including the frequency 
derivative in the light-curve folds, as each observation is 
sufficiently long that spin-down smears the pulsed signal (A 0 = 
0.11 and 0.32 cycles for the first and second epoch, respectively). 
The final period measurement used the iterated period derivative 
reported in Table 2. The uncertainties are estimated from a 
Monte Carlo simulation of the light curve using the method 
described by Gotthelf et al. (1999). The final power spectra 
are shown in Figure 2 (left). The derived physical parameters 
in Table 2 are from the relations E = Ajx 2 I P P~ 3 , where 
I = 10 45 g cm 2 , and the dipole spin-down relations z c = P /2P 
and B s = 3.2 x 10 l9 (PP) 1 ' 2 G. 

Figure 2 (right) shows the highly modulated pulse profile 
characterized by a relatively sharp peak and broad trough 
compared to a pure sinusoid signal. The pulsed fraction in 


Table 2 

Timing Parameters of PSR J1640— 4631 


Parameter 

Value 

R.A. (J2000.0) a 

16 h 40 m 43?52 

Decl. (J2000.0) a 

— 46°31'35."4 

Epoch (MJD TDB) 

56,466.0 

Frequency 15 ,/ 

4.843950957(40) Hz 

Frequency derivative 15 , / 

-2.290(10) x 10““ Hz s -1 

Period 15 , P 

0.206443048(33) s 

Period derivative 15 , P 

9.758(44) x 10“ 13 

Spin-down luminosity, E 

4.4 x 10 36 erg s -1 

Characteristic age, r c 

3350 yr 

Surface dipole magnetic field, B s 

1.4 x 10 13 G 


Notes. 

a Chandra position from Lemiere et al. (2009). 
b lcr uncertainties given in parentheses. 


the 3-25 keV band is f p > 48 zb 10% after allowing for the 
background level in the source aperture, estimated using data 
from a concentric annulus. This increases to f p ~ 82% after 
taking into account contamination from the PWN in the source 
aperture using the spectral results presented below. Here, we 
define the pulsed fraction as the ratio of the pulsed emission to 
the net source (background subtracted) flux. To determine the 
unpulsed level we take the average of the two lowest bins in the 
16-bin folded light curve of Figure 2. The pulse shape shows no 
variation with energy, within statistics. 

2.2. Spectral Analysis 

For spectral study we extracted photons from the observations 
listed in Table 1 using an elliptical region of diameter 3^4 x 
2^6 whose major-axis is oriented at position angle —40°, 
centered on the apparent enhanced PWN emission at (J2000) 
16 h 40 m 42. s 05, —46° 3 1/47 ."8, which is offset by 20" from the 
pulsar (see Figure 1). We estimate the background using nearby 
50" radius circular regions carefully chosen to account for stray 
light in the focal planes in two of the observations. Response 
matrices were generated for each spectral file using the NuSTAR 
analysis software. All eight spectra were combined, as were 
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Table 3 

X-Ray Spectrum of PSR J1640— 4631 and its Wind Nebula 


Parameter 

Chandra Only 

Chandra + NuSTAR 

A h (cm -2 ) 

(1.2 ±0.6) x 10 23 

(1.8 ±0.6) x 

10 23 

r PSR 

1 2 +0 - 9 
1 * Z -0.8 

1 3+0.9 
^ — 0.5 


F psr (2-10 keV) 

1.9t°i 2 4 x 10-13 

(1.8 ±0.4) x 

10 -13 

r pwn 

2 3 +1 - 2 
z * — 1.0 

r\ rt+ 0.7 

z - z -0.4 


Fpwn (2-10 keV) 

5 .4^ 2 ' 6 3 x 10“ 13 

(5.5 ±0.8) X 

10 -13 

Xv ( dof ) 

1.0(56) 

0.82 (79) 


Notes. Absorbed power-law model fits to the pulsar and the PWN spectra with 
their column densities linked. The simultaneous fit to NuSTAR and Chandra 
data is described in the text. The uncertainties are 90% confidence limits 
determined from the error ellipse extrema. The given fluxes are absorbed, in 
units of erg cm -2 s -1 . 


their matching response matrices, and fitted using the XSPEC 
software package (Arnaud 1996) with x 2 minimization. The 
fitting is limited to the 3-20 keV range due to low signal-to- 
noise at higher energies. In this range, the combined spectrum 
yields a total of 14,100 source plus background photons, and is 
grouped to obtain >5a signal-to-noise per spectral bin. 

The NuSTAR spectrum in the large aperture is dominated 
by PWN emission and is well fit by an absorbed power-law 
model as expected for non-thermal emission from the nebula. 
We used the tbabs absorption model, with Wilms et al. (2000) 
abundances and Verner et al. (1996) cross sections, and obtain 
a best-fit Ah = (1.7 zb 0.9) x 10 23 cm -2 and T = 1.9 ± 0.4. 
Errors are 90% confidence level (Ax 2 = 4.61 for two interesting 
parameters) determined from the error ellipse contour extrema. 
This provides an acceptable fit with a reduced x 2 = 0-87 for 
36 degrees of freedom (dof). The absorbed 2-10 keV flux 
is (8.O+20) x 10 -13 erg cm -2 s -1 (90% confidence). This 
is consistent with the values reported by Funk et al. (2007) 
using XMM-Newton data extracted from a 2^5 diameter aperture 
centered on the PWN. 

To further constrain the PWN emission and to isolation 
the pulsar contribution to the NuSTAR spectrum, we include 
the Chandra spectrum of the pulsar and PWN in the fitting 
process. We supplemented the Chandra data set (ObsID 7591) 
previously analyzed as part of the study of HESS J1640— 465 
by Lemiere et al. (2009) with data obtained on 2011 June 6 
(ObsID 12508). These data were acquired with the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) and 
reprocessed using the chandra_repro package of the Chandra 
Interactive Analysis of Observations (CIAO) software suite. 
Spectra and their response matrices from the two observations 
were produced using specextract. The pulsar and PWN 
spectra are grouped with a minimum of 10 and 40 counts 
per spectral bin, respectively, and each fitted with the absorbed 
power-law model, using a common column density. 

We used an extraction radius of 2 " for the pulsar; the 
background is negligible in this small region. For the PWN 
we used the NuSTAR elliptical region and similar background 
region. We combined spectra from the individual observations, 
and their response matrices, as above. Although the ACIS- 
I detector is sensitive in the 0.5-10 keV energy range, we 
restricted the spectral fits to the 3-7 keV range due to the high 
absorption and limited statistics. A total of 142 and 1 369 photons 
were fitted for the pulsar and PWN, respectively, during the 
47 ks exposure. The second column of Table 3 presents the 
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Figure 3. Chandra and NuSTAR spectra of PSR J1 640— 4631 (red) and its 
wind nebula (blue) fitted simultaneously with the absorbed power-law models 
presented in Table 3. The pulsar and PWN components of the NuSTAR spectrum 
sum to the total NuSTAR spectrum (black). The Chandra and NuSTAR fits are 
tied in the 2-10 keV band, but allow for a constant flux offset between missions. 
The lower panel shows residuals from the best fit. 

(A color version of this figure is available in the online journal.) 

resulting spectral fits for the pulsar and PWN spectra using the 
Chandra data alone. 13 

We fitted the Chandra and NuSTAR spectra together using 
two power laws as before, with their column densities tied. 
To constrain the pulsar and PWN contributions to the NuSTAR 
spectrum, we fixed the power-law indices to the Chandra models 
and tied their flux together in the overlapping 2-10 keV band. To 
allow for flux calibration differences between the two missions 
we introduce an overall normalization constant to the NuSTAR 
model, with best fit value of 1.11. The resulting spectrum is 
shown in Figure 3, and the parameters are reported in Table 3. 
The NuSTAR spectra of the pulsar and PWN are evidently 
successfully modeled, as the fitted values are in good agreement 
with the Chandra results for each component. To estimate 
the PWN contribution to the pulse profile shown in Figure 2, 
we repeated our spectral analysis using a 30' radius aperture; 
the PWN is found to account for ~53% of the background- 
subtracted source flux in pulse profile. 

3. SEARCH FOR GAMMA-RAY PUFSATIONS 

Although Slane et al. (2010) argue that the GeV emission 
from 1FGF J1 640.8—4634 originates from the PWN, the Fermi 
source is spatially unresolved and marginally consistent with 
the position of the pulsar. Therefore, we searched the Fermi 
FAT data at the position of PSR J 1640— 4631 for a pulsed 
signal. The X-ray timing parameters and their errors given 
in Table 2 allow a search for pulsations around the known 
values. We extract photon arrival times from 2008 August 4 
to 2013 October 2. These data were reprocessed using the 
“Pass 7” event reconstruction algorithm. We selected “source” 
photons for zenith angles < 100° and restricted the energy range 
to >500 MeV to minimize Earth limb and diffuse Galactic 
y-ray background contamination. The photon arrival times were 
corrected to the solar system barycenter using the JPF DE405 

13 The fluxes described in Lemiere et al. (2009) as unabsorbed are not 
consistent with the (absorbed) fluxes given in Table 3 or with those presented 
by Funk et al. (2007); the Lemiere et al. (2009) values are likely absorbed 
fluxes, mischaracterized. 
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ephemeris and the Chandra coordinates. For our nominal pulsar 
search we extracted photons from an energy-dependent radius 
enclosing 95% of the point spread function. 

Because the sparse X-ray observations do not provide a phase- 
connected ephemeris, unknown timing noise could limit the 
practical time span of a coherent pulsar search. Accordingly, we 
grouped the photons into intervals of 100 days and searched for 
a significant signal over a range of /, / centered on interpolated 
values at the test epoch using the X-ray timing parameters. This 
range corresponds to three times the la uncertainty on each off 
and /. We used the if-test statistic (de Jager & Biisching 2010) 
to allow for a complex profile with narrow features. This test 
selects the Fourier harmonic m ^ 20 that results in the most 
significant normalized power. 

None of the resulting searches yielded a significant signal 
above the noise. We also summed the power from these 
searches for each /, / pair, to increase the signal sensitivity. 
No significant signal was apparent. To explore a range of 
instrumental parameters that might be more sensitive, we 
repeated our search for a low-energy cutoff of 100, 300, and 
500 MeV, and time interval of 30, 200, and 300 days. Despite 
the expanded range of parameter space, no signal stronger than 
the expected noise is detected in each search or in the summed 
results. 

4. MODELING THE BROADBAND SPECTRUM 

4.1. Age, Energetics, and Distance 

The discovery of PSR J1640— 4631 supports the conjec- 
ture that CXOU J164043. 5— 463135, HESS J1640-465, and 
1FHL J1640.5— 4634 are all manifestations of a middle-aged 
pulsar/PWN system. For the measured distance of 12 kpc to 
G338.3— 0.0 (see below), the 0.2-10 TeV luminosity of HESS 
J1640— 465 is 6% of the pulsar’s present E , while the ratio 
E(0.2-10TeV)/F PWN (2-10keV) ~ 13, indicating that, in a lep- 
tonic model where the TeV emission results from up scattering 
of ambient photons by relativistic electrons, inverse Compton 
losses now dominate over synchrotron emission. 

A result that was perhaps unexpected is the young char- 
acteristic age of the pulsar. The pulsar’s characteristic age, 
t c = P/2P = 3350 yr, is an approximate measure that ap- 
plies when a simple vacuum dipole spin-down model is as- 
sumed, and when the spin period is significantly larger than 
that at birth. Although pulsars for which the model has been 
tested do not rotate like simple vacuum dipoles (Kaspi et al. 
2001), it is notable that, excluding magnetars, there are only 
seven pulsars with r c < 3000 yr in the Australia Telescope Na- 
tional Facility (ATNF) catalog 14 (Manchester et al. 2005), fol- 
lowed by PSR J1640— 4631 with z c = 3350 yr. This is younger 
than most authors have assumed in modeling the evolution of 
HESS J1640— 465/G338.3— 0.0 (T = 20 kyr, Funk et al. 2007; 
T = 15 kyr, Lemiere et al. 2009; T = 10 kyr, Slane et al. 
2010). Fang & Zhang (2010) presented two models similar to 
that of Slane et al. (2010), one with E = lx 10 38 erg s _1 and 
T = 4500 yr, and the other with E = 1.65 x 10 37 erg s _1 and 
T = 8200 yr. While these ages are closer to the likely true age, 
the actual spin-down power of PSR J1 640—463 1 is considerably 
smaller than the assumed values. 

Prior models for G338.3— 0.0 estimate the age of the SNR 
from its radius assuming that it is in the Sedov phase with 
r s = (2.02 S/p) l/5 t 2/5 (Spitzer 1978). These estimates are 


14 http://www.atnf.csiro.au/research/pulsar/psrcat/expert.html 


imprecise due to the unknown energy S of the explosion and 
the ambient mass density p, and the strong dependence of 
age t on radius, which magnifies any uncertainty in distance. 
Also, the pulsar is offset from the center of the SNR in a 
direction that, given the morphology of the PWN, cannot be 
explained by high kick velocity. Lemiere et al. (2009) conclude 
that the geometrical center of the remnant is not the explosion 
site. Instead, the structure of the remnant may be affected by 
complicated interactions with the local ISM, which could affect 
the estimation of its age. 

There is no near/far distance ambiguity for G338.3— 0.0 
because 21 cm absorption is seen up to the maximum negative 
value at the tangent point in both it and adjacent Hu regions. 
Lemiere et al. (2009) concluded that the distance to G338.3— 0.0 
is in the range 12-13.5 kpc, while Kothes & Dougherty (2007) 
had already found a value of ~ 1 1—12 kpc from the same 21 cm 
data, but assuming a Galactic center distance of 7.6 kpc instead 
of 8.5 kpc. We adopt d = 12 kpc here. While the distance 
appears to be constrained well, there is no direct measurement 
of the ambient density; the above referenced models assume 
values of tzism between 0.1 cm -3 and 10 cm -3 . Castelletti et al. 
(2011) estimated an electron density of 100-165 cm -3 for the 
western and northern part of the SNR shell by assuming that the 
low-energy turnover of its radio spectrum is due to local free- 
free absorption. They also searched the NANTEN CO survey 
data for molecular gas that could be target material for hadronic 
production of y-rays from protons accelerated in the SNR shell, 
but did not find any that could be associated with G338.3— 0.0. 

4.2. A Hadronic Model 

HESS Collaboration (2014) present a model in which the 
TeV emission is produced by hadronic interactions in the north 
and west part of the SNR shell. It is motivated largely by the 
continuity of the Slane et al. (2010) Fermi spectrum with the 
HESS spectrum, which they argue is difficult to fit in a PWN 
model, and also by the overlap of the TeV source with the side 
of the SNR shell that is adjacent to H ii regions. They require a 
high ambient density of 150 cm -3 to reproduce the GeV/TeV 
spectrum, but a small density to account for the size of the 
SNR shell. These are achieved by assuming that the explosion 
occurred inside a wind-blown bubble of the progenitor star, of 
average density 0.1 cm -3 , which is then consistent with their 
assumed age of 2500 yr. As they remark, the model requires 
a large fraction of the SN energy to be channeled into cosmic 
rays, especially since only one side of the SNR contributes to 
this process. 

As we argue below, the Slane et al. (2010) Fermi spectrum 
used by HESS Collaboration (2014) is not likely to be correct, 
and there may not even be a detectable source in Fermi at 
energies <10 GeV. This is why we use different Fermi results in 
our own model, which contributes to the diverging conclusions 
of the two papers. 

4.3. Leptonic PWN Models 

Slane et al. (2010) modeled the y - ray emission from HESS 
and the associated broad-band spectrum with an evolving, one- 
zone PWN model, where the y-rays are ambient photons inverse 
Compton scattered by the relativistic electrons that also produce 
the synchrotron emission. Crucial to these models are the age 
and spin-down power of the pulsar. Slane et al. (2010) assumed 
a remnant age of 10 kyr and a pulsar with spin-down power of 
E = 4 x 10 36 erg s _1 . The latter came from the E/L x relation 
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Table 4 

Properties of HESS J1640-465 


Gotthelf et al. 


Parameter 

Observed 

Modeled 

Reference 

SNR radius 

4:5 ± o:5 

4.'51 

Shaver & Goss (1970) 

S v [660 MHz] (mJy) 

^690 

337 

Giacani et al. (2008) 

Fx [2-25 keV] (erg cm -2 s -1 ) 

(1.68 + 0.4) x 10“ 12 

1.0 x 10“ 12 

This work 

T x 

9 9+0.7 
Z,Z -0.4 

2.4 

“ 

F y [10-500 GeV] (erg cm“ 2 s -1 ) 

(3.15 ± 1.00) x 10-“ 

1.6 x 10- 11 

Ackermann et al. (2013) 

u 

1.92 + 0.24 

1.96 



Note. The quoted 2-10 keV flux is corrected for interstellar absorption. 


of Possenti et al. (2002) and, given the observed scatter, it was 
fortuitously a very accurate prediction. It was necessary to add a 
Maxwellian distribution of electrons around 0. 1 TeV to a power- 
law tail to account for the strong Fermi GeV y - rays that they 
found from this source relative to its TeV flux. 

The pulsar age assumed by Slane et al. (2010) was likely too 
large. The true age of a pulsar is given by 


T = 


P 

(n - 1 )P 




where n = f f / / 2 is the braking index, assumed by Slane 
et al. (2010) to be 3, consistent with a vacuum dipole, Pq is the 
spin period at birth, and P is the currently measured period. For 
n — 3 the pulsar’s true age is less than its characteristic age, 
and refined models will likely favor a younger system and/or a 
smaller braking index. 

Assuming again that the PWN of PSR J1640— 4631 is pri- 
marily responsible for the y - ray emission, we fit an evolutionary 
model of a PWN inside an SNR to the broad-band spectrum. 
This can constrain the properties of the central neutron star, the 
pulsar wind, progenitor supernova, and the surrounding ISM. 
We include the PWN emission in the X-ray band and the radio 
upper-limit reported in the literature (see Table 4). However, 
in examining more recent Fermi publications, we note that the 
source 2FGL J1640.5— 4633 (Nolan et al. 2012) is fainter than 
the flux (based on ~ 1 yr of data) that was extracted by Slane 
et al. (2010). More recent analyses of this Fermi source have 
been restricted to energies >10 GeV (Acero et al. 2013; Ack- 
ermann et al. 2013), and we use only the latter in our revised 
modeling, in particular the spectrum of 1FHL J1640.5— 4634 
(Ackermann et al. 2013) derived from three years of data. In- 
deed, it is not obvious from inspecting Fermi counts maps that 
there is a significant source at energies <10 GeV, which is why 
we do not attempt to model this part of the Fermi band. This is 
an important difference from HESS Collaboration (2014), who 
accept and model the Fermi results of Slane et al. (2010). 

Our model is based on the work of Gelfand et al. (2009), 
modified to include background photon fields in addition to the 
cosmic microwave background (CMB), and a broken power-law 
spectrum of particles injected into the PWN at the termination 
shock, as favored by recent studies of these objects (e.g., 
Bucciantini et al. 201 1). We assumed two additional background 
photon fields: one with temperature 7j = 15K and energy 
density u\ — 4ucmb , where uqmb = 4.17 x 10 _13 erg cm -3 is 
the energy density of the CMB, and the other with T 2 = 5000 K 
and U 2 = 1.15i/cmb, the same as used by Slane et al. (2010). 
We fixed the distance at 12 kpc. To obtain the best- fit values of 
the free parameters and their errors, which are listed in Table 5, 
we employed a Markoff Chain Monte Carlo algorithm similar 


Table 5 

PWN Evolutionary Model Results for HESS J1640— 465 


Parameter 

Value 

SN explosion energy, log (£/10 51 erg) 

-0.78 

Ejecta mass, log (M e j/M 0 ) 

-0.34 

ISM density, log (nisM/cm -3 ) 

-1.53 

Braking index, n 

+1.87 

Spin-down timescale, log (to /yr) 

+2.89 

Wind magnetization, log (jib) 

-1.76 

Low-energy particle index, p\ 

+0.28 

High-energy particle index, p 2 

+2.59 

Minimum energy, log (E m \ n /GeV) 

-3.86 

Break energy, log (Efr/TeV) 

-0.71 

Maximum energy, log (£ max /PeV) 

+0.12 


to the one used to fit the properties of the PWNe in Kes 75 
(Gelfand et al. 2014) and G54. 1+0.3 (J. D. Gelfand et al. 2014, 
in preparation). 

For a constant braking index, the spin-down power of the 
pulsar evolves as 


E(t) = E 0 



-(^i) 


where Eo is the initial spin-down power and To is the initial 
spin-down timescale. The present P and P are known, which 
fixes the values of Eq and the present age T corresponding to 
each trial value of n and To. 

As shown in Figure 4, the model is able to reproduce the 
observed broad-band spectrum with y 2 = 24.7 for 15 degrees 
of freedom, using the best-fit parameters given in Table 5. In 
general, the values of these parameters are similar to what 
has been inferred for other systems using similar models (e.g., 
Bucciantini et al. 2011). However, our best- fit model requires 
a small braking index n ~ 1.9 and a short initial spin-down 
timescale To ~ 800 yr (see Figure 5). This implies that the pulsar 
was bom spinning rapidly, with Po ~ 15 ms and a high initial 
spin-down power Eq ~ 10 40 erg s -1 . Its present age would then 
be 6800 yr. We note that, in this case, the model is not entirely 
self-consistent, as the initial rotational energy of the pulsar, 
9 x 10 49 erg, is a significant fraction of the fitted SN explosion 
energy, and could affect the dynamics of the explosion. 

5. DISCUSSION 

A PWN model for HESS J1 640—465 is easier to accom- 
modate now because we have deprecated the original Fermi 
spectmm that was used in all previous model fitting. Even 
though that spectrum appeared to be continuous with the new 
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Figure 4. Broadband PWN model fit to the spectrum of HESS J1 640— 465. From 
left to right, in red, the radio upper-limit (Giacani et al. 2008), the Chandra X-ray 
spectrum of the PWN (Lemiere et al. 2009), the Fermi GeV spectra (Ackermann 
et al. 2013), the HESS TeV data points (HESS Collaboration 2014). The blue 
lines represent the error range on the Chandra and Fermi spectral fits. The black 
line is the spectral energy distribution for the best-fit model parameters given in 
Table 5. 

(A color version of this figure is available in the online journal.) 

HESS points, fitting a single power law, it is now evident that 
additional Fermi data above 10 GeV are continuous with the 
HESS data, but with spectral curvature to lower energy that is 
characteristic of lepton cooling. As to the spatial distribution 
of TeV photons, while they overlap with part of the SNR shell 
of G3 3 8. 3— 0.0, they also overlap with the X-ray PWN, so it 
is entirely possible that both the PWN and the shell contribute 
to the y-ray emission. Furthermore, it cannot be ruled out that 
PSR J 1640— 4631 emits some pulsed y-rays around 1 GeV that 
we have not yet been able to discover. 

The modeled SN explosion is energetic enough for 
G3 3 8. 3— 0.0 to have reached its present radius of 14 pc in 
6800 yr. The fitted ambient density is ~0.03 cm -3 , similar to the 
value assumed by HESS Collaboration (2014), so this parameter 
is not a discriminant between the models. The rapid spin-down 
of PSR J 1 640—463 1 helps to explain why its PWN can be such a 
luminous y- ray source. The electrons currently emitting inverse 
Compton scattered TeV photons were injected when the pulsar’s 
E was much higher, while the PWN has expanded, reducing its 
magnetic field strength and limiting its synchrotron losses. A 
similar scenario may explain why the young magnetars SGR 
1806-20 (Rowell et al. 2012) and CXOU J171405.7-381031 
in the SNR CTB 37B (Halpern & Gotthelf 2010) may power 
TeV sources even though their present spin-down luminosities 
are small. 

If n < 3 for PSR J1 640—463 1 , as it is for all pulsars in which 
it has been measured, one reason could be that its magnetic 
dipole field strength or the inclination angle a of its axis is 
increasing in time such that ( B s sin a) 2 oc P 3 ~ n (Livingstone 
et al. 2007). If n = 2, for example, P will remain constant and 
PSR J 1640— 4631 will join the group of magnetars when its 
period is ~7 s and B s = 8 x 10 13 G. (However, aside from the 
idealized n = 3 case, there is no reason why n should be constant 
over the lifetime of a pulsar.) The smallest measured pulsar 
braking index that is not affected by glitching isn = 0.9±0.2 
for PSR J1734— 3333 (Espinoza et al. 2011). The 1.17 s period 
of PSR J1734— 3333, and its inferred dipole magnetic field of 
~5 x 10 13 G, are already close to those of magnetars. 


Figure 5. Trial values of the braking index n and spin-down timescale to of PSR 
J1640— 4631. The grayscale indicates the distribution of trials with x 2 > 27.7, 
while the colors indicate the x 2 of the particular trials with x 2 ^27.7. 

(A color version of this figure is available in the online journal.) 

The complex of H n regions and SNRs including G338.3— 0.0 
is at the intersection of the far end of the Galactic bar with the 
Norma spiral arm. Davies et al. (2012) detected a young, massive 
star cluster in the near-infrared at the center of the H ii region 
G338.4+0.1, and displaced by 8' from HESS J1640-465, or a 
projected separation of 22 pc at the distance they inferred of 
1 1 zb 2 kpc. The most massive star in the cluster is a WR star 
with an estimated initial mass of 62 M 0 and an age of 3.7 Myr. 
The optical extinction to the cluster is compatible with the 
X-ray measured column density to CXOU J164043. 5— 463135. 
Optical and infrared emission from the Hii regions and star 
cluster could provide significant target photons for inverse 
Compton scattered TeV emission. Davies et al. (2012) suggested 
that the progenitor of the putative neutron star could have been 
born in, and dynamically ejected from the cluster during its 
formation, in which case its mass must have been at least as 
large as the most massive star presently in the cluster. But they 
stopped short of making this claim because they could not prove 
the birth of the pulsar progenitor in the cluster, as opposed to in 
a nearby site of independently seeded star formation. 

Nevertheless, it is possible to speculate that the large magnetic 
field strength of PSR J1640— 4631, and the short initial spin 
period needed to explain its powerful TeV nebula, may both 
result from a massive progenitor such as a WR star. Duncan & 
Thompson (1992) proposed that millisecond initial spin periods 
are needed for a dynamo to generate magnetar strength 5-fields. 
There is evidence that the most massive stars (>40 M 0 ) are 
the progenitors of magnetars (Figer et al. 2005; Gaensler et al. 
2005; Muno et al. 2006), although there are exceptions (Davies 
et al. 2009). If PSR J1640— 4631 is found to have these birth 
properties, it could be evidence of a physical link between 
magnetars and other high 5-field pulsars. 

6. CONCLUSIONS AND FUTURE WORK 

NuSTAR is a sensitive instrument for the discovery of young 
pulsars in distant parts of the Galaxy that are obscured by 
large ISM column densities. The detection of PSR J1 640— 4631 
in G3 3 8. 3— 0.0 provides long-awaited evidence that a PWN 
powers the TeV source HESS J1640— 465, and its properties 
test evolutionary models that fit the multiwavelength spectrum. 
The spin-down power of PSR J1640— 4631, 4.4 x 10 36 erg s -1 , 
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is typical of the range of middle-aged pulsars powering TeV 
nebulae, but its characteristic age of 3350 yr is younger than 
the modeled age of the system, which is ~6800 yr, with a 
corresponding braking index of n ~ 2. Its observed rapid spin- 
down and inferred short initial spin period of Pq ~ 15 ms 
make it possible for PSR J1640— 4631 to power the most 
luminous known TeV source in the Galaxy. Nevertheless, a more 
sensitive determination of the y-ray spectrum below 10 GeV, 
and a map of the spatial distribution of the TeV photons, are 
needed to discriminate among models that locate the emission 
in the PWN versus the SNR shell. It is possible that two (or 
more) mechanisms may contribute to the extraordinary y - ray 
luminosity of HESS J1640— 465. 

If radio pulsations are detected and monitored, or with the 
dedicated program of timing PSR J1 640—463 1 with NuSTAR in 
progress, we can hope to measure its braking index, and thereby 
better estimate the actual age of G3 3 8. 3— 0.0 and the initial spin 
period of the pulsar. This can provide confirmation of the model 
presented in this work. 

This work made use of data from the NuSTAR mission, a 
project led by the California Institute of Technology, managed 
by the Jet Propulsion Laboratory, and funded by the National 
Aeronautics and Space Administration. This work has also made 
use of archival data from the Chandra X-ray Observatory. Par- 
tial support for E.V.G. and J.A.T were provided by NASA 
through Chandra Award Numbers S AO G02- 1 3097X and GO 1 - 
12068 A, respectively, issued by the Chandra X-ray Observa- 
tory Center, which is operated by the Smithsonian Astrophys- 
ical Observatory for and on behalf of NASA under contract 
NAS8-03060. E.V.G acknowledges support by Fermi Guest In- 
vestigator Grant NNX12A089G. Radio contours were obtained 
from Molonglo Observatory Synthesis Telescope (MOST) data 
provided by the University of Sydney with support from the 
Australian Research Council and the Science Foundation for 
Physics within The University of Sydney. 
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